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Abstract We consider the Nielsen identities for the two-point functions of fuU QCD 
and QED in the class of Lorentz gauges. For pedagogical reasons the identities are first 
derived in QED to demonstrate the gauge independence of the photon self-energy, and 
of the electron mass shell. In QCD we derive the general identity and hence the identi- 
ties for the quark, gluon and ghost propagators. The explicit contributions to the gluon 
and ghost identities are calculated to one-loop order, and then we show that the quark 
identity requires that in on-shell schemes the quark mass renormalisation must be gauge 
independent. Furthermore, we obtain formal solutions for the gluon self-energy and ghost 
propagator in terms of the gauge dependence of other, independent Green functions. 

PACS No.: 11.30.Ly 12.38.Aw 12.20.Ds 



1. Introduction 

The Nielsen identities'^] deserve to be better known. They follow from a modification of 
Becchi-Rouet-Stora-Tyutin (BRST) symmetry in a similar way to how the usual Ward 
identities are obtained from BRST. The identities are, however, of a quite different nature 
to the Ward identities. This may be illustrated by (very) schematically writing them as 
-§fGi = 6*2^^35 where ^ is the gauge parameter and the Gi are Green's functions. This 
then, for example, clearly tells us that if the right hand side vanishes, then Gi must be 
gauge independent. The aim of this paper is to investigate the consequences of the Nielsen 
identities for the two-point functions of QED and QCD. 

We will not follow here the original derivation of Nielsen, but rather use an alternative 
approach due to Piguet and Sibold'^l. The identities have been previously employed ["^^ to 
study the gauge invariance of the Higgs mass and the role of gauge symmetry in the 
effective potential for models of spontaneous symmetry breaking. 

It will become evident in what follows that the Nielsen identities are particularly 
useful for investigations of on-mass shell Green's functions and on-shell renormalisation 
constants. The on-shell renormalisation scheme is widely used in QED, the electroweak 
theory and QCD with heavy quarks and it is of the greatest importance for studies of the 
S-matrix. 

For completeness and clarity we begin by discussing the case of QED in Section 2. We 
show in particular that the pole mass of the electron is gauge independent and that the 
photon self-energy can be simply shown to be gauge parameter independent. We stress 
here the interplay between the usual (BRST) Ward identities and the Nielsen identities. 

In Section 3 we present the Nielsen identities for QCD. We explicitly calculate the 
one-loop contributions in the identity for the gluon propagator, and then demonstrate that 
in on-shell renormalisation schemes the quark mass renormalisation must also be gauge 
independent to all orders. The Nielsen identity for the ghost propagator is constructed 
and its explicit content is determined to one-loop order. Formal solutions for the gluon 



self-energy and ghost propagator appear from this analysis which deserve further study. 

Section 4 provides a perturbative analysis of the Nielsen identities illustrating the 
gauge parameter independence of the mass renormalisation in the mass-shell scheme. Some 
conclusions are presented in Section 5. 

2. The QED Identities 

The standard path integral formulation of QED in covariant gauges is based upon the 
following Lagrangian 

-Cqed = -\f^ + i^iip- rn)ij + ^B'' + Bd-A- cDc, (2.1) 

where ^ is the gauge parameter, F is the field strength and B is an auxiliary field ^'^^ 

Although ghosts decouple in QED, one often includes them so as to demonstrate that the 

Lagrangian is invariant under the following BRST transformations 

5A^ — ed^c , di/j = +iec'i/j , 5c = , 

(2.2) 
SB = , difj = —ieci/j , dc = eB . 

The Ward identities may then be obtained by exploiting this invariance in the usual 

fashion [^]. 

Following [2], the Nielsen identities are obtained by making the following addition to 
the Lagrangian 

-Cqed -^ -Cqed + —cB , (2.3) 

where x is a global Grassmannian variable, x^ = 0. (In the following care must be taken 
with the minus signs that are needed under the interchange of Grassmannian variables!) It 
is clear upon a little reflection that the addition of this term cannot change the dynamics 
of the theory. To see this most easily consider that the generating functional can be 
expanded in x and, as a result of its Grassmannian nature, only two terms survive: x^ 
and x^. The first term gives us the usual dynamics of QED and the second must vanish 
by virtue of ghost number when we integrate over the ghost fields. This shows that we 
have not changed any physics. In what follows we will employ the Lagrangian (2.3). 
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This modified Lagrangian is invariant under the foUowing extended set of BRST trans- 
formations 

(5+A^ = ea^c, S+c = eB, 

6+B = 0, 6+c = 0, 

(2.4) 
S+if; = -lec^ , 5^i = ex , 

d'^i/j = +iec'i/j , d~^x = • 

where e is a Grassmann quantity. The F^ and Dirac Lagrangians are, as usual, separately 
invariant and the remaining terms are readily seen to be invariant. The crucial point to 
note is that the gauge parameter is now transformed. 

To exploit this invariance and so derive the Nielsen identities we now consider the 
following generating functional 

Z = [dn] exp < i jct^x £qed + J^^A^+J-^,^ + ^J^ + BJb 

J ^ J ^ (2.5) 

+ K^{-icip) + icipK^ 

The various sources denoted by J with a subscript are standard ones. We have not included 
sources for the ghost fields, for example, since we will not consider their Green's functions, 
for QED they are anyway trivial. The purpose of the additional, rather exotic looking, 
sources denoted by iiT's will become apparent in a moment. Note that we may rewrite this 
part of the Lagrangian as 

To study the gauge dependence of the electron and photon propagators, we now 
introduce the generating functional of proper Green functions. 

M^( J^, J^,J^, Jb, K^, K^, X, - jd'^x J^A^ + JbB + J^^ + ^J^ 
As a consequence of invariance under (2.4) we have 



(2.7) 



where we have used S^B = 5+c = S~^x = 0. From (2.6) we have S'^^i/j = e-~— etc. We may 
therefore rewrite (2.8) as 

^ 5V 5V 5V 5V 5V ^5V dV , , 

oA^ 6K^ Qi\) oK^ St/j be oi 

Note that imphcit coordinate space integrations are understood in the above! 

We are now in a position to obtain the Nielsen identities for QED. We differentiate 
(2.9) with respect to x and then set x fo zero. 

H H = + S 

SxSK^ 5ip SK^ 5x5ip 5x5c 

In the special case when no further functional derivatives with respect to ghost fields will 
be applied to (2.10a) ghost number conservation implies a simplification. The so simplified 
result is then 

dV ^5^V d^V dV d'^V dV , ,, 

where we have used the result d^c = d^jj: = for the one-point function. From these 
central results we can generate the QED Nielsen identities for the two-point functions. 

2.1 The Mass of the Electron 

The first application studied here is the electron mass. Although the the fermion field 
ip is not BRST-invariant and may not be naively identified with the electront^'^l, we will 
now show that its pole mass is gauge independent and may so be given a physical meaning. 
The inverse fermion propagator is given by iS~^{y—x) = ^,/.F,, . . Differentiating (2.10b) 
with respect to V'(y) ^^^ i/j{x) we obtain 

f ,,, ,,,> , ,7, ,,, + ,,, ,,7 ,7, ,,r. . ' (2-11) 



d^ 6ijj{y)6ijj{x) 6ijj{y)6K^6xSipix)Sijj 6ip{y)6ijj 6ip{x)6K^6x 
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and all other terms will vanish, either by fermion conservation or by their relation to 
one-point functions. 

To complete the analysis of the Fermion propagator (2.11) must be transformed into 
momentum space. Defining 

d^r f dS dH 

(2.12) 



dxSK^H^ijiy) J (27r)4(27r)4 

^^r f d^q d^e 






5xdK^{w)5^{y) J (27r)4(27r)4 
and recalling that (2.11) contains implicit w and z integrations, it is found that 

d 



d^ 



S-\p) = S-\p){F{p,-p,Q) + F{-p,p,Q)] 



(2.13) 



This result is the Nielsen identity for the inverse fermion propagator. In particular, since 
the right hand side vanishes at the mass shell, and since F(— p, p, 0) has no single particle 
poles, we have 

0, (2.14) 



dS-^{p) 



di 



p2=M'2 



where M is the pole mass. 

From this result it is easy to see that M is gauge independent. Since we may, in a 
covariant gauge, generally write the inverse propagator as 



S-\p) ^ Aip^)i>- Bip' 



(2.15) 



then M is defined by 

A{M^)M = B{M^) . (2.16) 

If we differentiate this equation with respect to ^ and compare it with (2.14) rewritten in 
the following way 



f dAjp'') dBjp''] 



^^ dA{M^) dBjM^) ^ 



p'2=M'2 



di 



we see that we obtain the desired result 



A + 



dA dB 



dM dM 



dM 



dB{M^) dA{M^] 



dC 



dC 



di 



dM 

^ -r- =0 



(2.17) 



di 



(2.18) 



so that the pole mass is gauge parameter independent (note that the bracketed quantity 
on the left hand side of (2.18) is non-zero). In Section 4 this point will be studied more 
carefully in a perturbative fashion. 

2.2 The Photon Propagator 

The inverse photon propagator may be studied in a similar fashion to the above. By 
functional differentiation of (2.10b) with respect to Aiy{x) and Ax{y), we find 

^ ^!E = (2 19) 

dCSA,ix)6Ax{y) ' ^ ■ ' 

where many terms that must, from considerations of fermion and ghost number or relation 
to one-point functions, vanish have been neglected. At this stage it is important to recog- 
nize that the auxiliary field B is independent of A^. To make a direct connection between 
(2.19) and the photon vacuum polarization, it is necessary to consider some aspects of the 
auxiliary field formalism [^'^1. 

To formulate perturbation theory it is necessary to consider the free field case, corre- 
sponding to the quadratic part of the Lagrangian (2.1). The mixing between B and d-Ain 
the Lagrangian must be diagonalized after functional integration, leading to the following 
free Geld bosonic propagators. 

f d^x e^P-^(0|T {B{x)B{0)) \0) = (2.20a) 

d^x e*P-"(0|T (B(x)AJO)) \0) = ^ (2.206) 



d'^x a'P-"" {0\T {A^{x)A^{0)) \0) = i 



-^ + (1-0^ 

pZ p'i 



(2.20c) 



BRS symmetry implies that (2.20a) and (2.20b) are true to all orders in perturbation 
theory. 

0= ^ (0|T {c{x)B{y)) \0) = {0\T {B{x)B{y)) \0) (2.21a) 

= ^(O|T(c(x)A,(y))|O) = (O|T(S(a:)A,(y))|O) + (O|T(V(x)c(y))|O)(2.216) 
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By recognizing that the ghosts are decoupled, (2.21b) then imphes that (2.20b) is vaUd to 
all orders in perturbation theory. This then implies that radiative corrections to the tree- 
level result (2.20b) are absent, so the photon vacuum polarization is necessarily transverse 
as illustrated in Figure 1. 



= p^ll^x 
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u\ 



+ (1-0 



p'^p^ 



p^ p^ 

Hence we have the usual result for the full photon propagator. 



(2.22) 



d^'x e^P-(-^) {0\T {A^{x)AM) \0) = i 
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flU 



pz 



+ 



p4 



1 



1 + n(p2 



-ti 



p4 



(2.23) 



The implications of the results for the Green functions (2.20a), (2.20b) and (2.23) for 
the generating functional F can be understood by consideration of the following functional 
identities. 



5B{x) 
5B{x) 

my) 

5A^{x) 







SAf,{x) 
6B{y) 

5{x-y) 
S{x-y)Si^ 



(2.24) 



SA^iy) 

In terms of generating functionals these expressions become 

6^r 6^W 6^r 



5{x - y)6i^ = / d'^z 



+ 



d^W 







d{x-y) 



d^z 



d^z 
d^z 



5Ax{z)5A^{y) 5Jx{z)5J^{x) 5B{z)5A,{y) 5Jb{z)5J^{x) 

dAx{z)dB{y) 6Jx{z)6J^{x) ^ 6B{z)6B{y) dJB{z)6J^{x) 

+ 



(2.25) 



6Ax{z)6A^iy) 6Jx{z)6Jb{x) 6B{z)6A^{y) 6Jb{z)6Jb{x) 

+ 



5Ax{z)5B{y) 5Jx{z)5Jb{x) 5B{z)5B{y) 5Jb{z)5Jb{x) 
The functional derivatives of W are known from the Green functions discussed above, and 
hence (2.25) implies that the quadratic pieces of F simply diagonalize the various A^, B 
two-point Green functions in momentum space. 

5^T 



d xe 



ip-x 



6A^{x)6A,{0) 



p'[-9^.u+P^.pA (l + n(p2)) 



(2.26a) 



X2p 

The transversality of (2.26a) is easily seen to be consistent with the fundamental 
identity (2.9). f Upon functional differentiation of (2.9) with respect to A^ and c, setting 
^ = 0, and then Fourier transforming, we see that 

p'' I d^xe'P-''—-^^-—- = Q (2.27) 

consistent with the analysis leading to (2.26a). 

Finally, substituting the result of (2.26a) into the identity (2.19) we find that the 
polarization 11 (p^) must be independent of the gauge parameter!^] 

—^=Q. (2.24) 

This result shows the power of the Nielsen identities. 

Since the ghosts decouple this concludes our survey of the Nielsen identities for the 
two-point functions of QED. 

3. The QCD Identities 

The reader who has closely followed the above will have no difficulty in finding the analo- 
gous identities for QCDt^L We give the basic steps and so define our notation. (Note that 
colour indices are implicit.) The modified QCD Lagrangian is 

>Cqcd = -\f^ + ^{ip- m)^ + 1^2 + sa • A - cd^D^ c + |cS (3.1) 

which is invariant under the following augmented BRST transformations 



5+S = 0, 5+c^-]-e[c,c] 

d'^i/j = —iecifj , d^^ = ex , 

5"'~'i/' = -t-zecf/j , d~^x = • 



(3.2) 



t We are grateful to the referee for bringing this point to our attention. 
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As in the QED case, the extension of BRST includes a transformation of the gauge pa- 
rameter. 

The generating functional Z with sources (including ghosts) is 
Z = / \d[}\ exp < i Id'^x Cqcd + J^iA^ + JbB + J^^j + ^J^ + JcC + cJc 

+ K^{D^c) -]^K,[c,c] + K^i-icij) + icijK^ \ . 
and the generating functional F for the proper Green functions is 

r(A^, -0, •0, s, c, c, X, ^, K^, K^, K^) =w{j^, Jb, J^, J^, Jc, Jc, k^, k^, i?^, X, 

- /d^x J^A'' + JbB + J^i; + ^ J^ + J^c + cj 



(3.3) 



(3.4) 



From this last result it is simple to obtain the nonabelian version of (2.9). 



After differentiation with respect to x ^-nd then setting % = we find 

aF d'^V dV _ dV d'^V d'^V 5V 

6r 6^r 6^r 6r 6r 6^r ^6^r , ,, 

+ 77T^T7 + .r^ . .T +Bj-j^^ (3.56) 



6K^ 6i)6x Sx6K^ 6%l) 6K^ 6x5'il^ ^X^c 
6'^V 6V 6T 6'^V 



6x6Kc 6c 6Kc 6x6c 

which is the QCD equivalent of (2.10a). As before, further functional differentiation with 
respect to fundamental non-ghost fields can be obtained from the following simpler result 
using ghost number conservation in (3.5b). 

aF 6'^V 6V 6'^V 6V 6'^V 6V 6^V 

" a^ + 6x6K^ 6A^ + '6^^J^ + '6^^M ^ ^^ ^ ' ""^ 

Although expressions (3.5c) and (2.10b) are superficially very similar, the distinction 
between abelian and non-abelian theories shows up in two important ways. First, there is 
now the current K^ which couples ghosts and gauge fields through the covariant derivative 
and secondly the ghosts are no longer decoupled. 
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3.1 Nielsen Identity for the Gluon Propagator 

The distinction between non-abelian and abelian theories is most evident when considering 
the Nielsen identity for the gauge field propagators. To our knowledge, these identities have 
not been previously studied in an explicit calculation. Proceeding as in Section 2.2, one 
finds 

+ z/ ^^ A . {6.oa) 



d^ 6A,{x)6Ax{y) SK^SxSA,{x) SA^6Ax{y) 
Defining the Green functions, 

we have the following result after Fourier transforming (3.5a). 

— r,A(p) = rA^(p)F^.(-p,p,0) +z/ ^ A (3.6rf) 

As in the photon analysis, it is necessary to relate the gluon vacuum polarization to 
^ ^lviJ>) before studying (3.6d). The analysis is only slightly more complicated than in the 
abelian case. The free field results (2.20) clearly are identical, but the BRS symmetry is 
somewhat diflFerent. 

= ^(0|T(c(x)S(y)) \0) = {0\T {B{x)B{y)) \0) (3.7a) 

= ^(O|T(c(x)A^(y))|O) = (O|T(S(x)A,(y))|O)-(O|T(c(x)I),c(y))|O)(3.76) 

The relevant Green function in (3.7b) is 

d^xe'P-^OYTicmD^cix)) \0) = ^ (3.7c) 

as can be verified by contracting both sides with p^ then applying the ghost equation of 
motion and the canonical commutation relations. 

Thus we see from (3.7) that the abelian results of (2.20a) and (2.20b) apply to all 
orders in QCD. The same argument as outlined in Section 2.2 then implies that the gluon 
vacuum polarization must be transverse, and so (2.26a) is also valid in QCD. 
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In contrast to the case of QED, (3.6d) cannot be solved exactly since the Green 
function F^j^(p, — p, 0) must be determined perturbatively. However, we can explicitly 
display the perturbative contributions in this identity, specifically to one-loop order. Since 
the vacuum polarization n(p^) is well known (see, e.g., [8]) we only have to calculate 

When calculating F^i,{p, —p, 0) it is important to recognize that the mixing between 
the fields B and A^ implies that we do not simply have a truncated Green function. This 
can be seen in the following way: 

S^T f ,4 S^W 6Jx(s) S^W SJb(s) 

~ a s ^^^ , , , — , ^ , , — - — -— + 



5Af,{x)6K^{y)6x J 6K^{y)6xSJx{s) SA^{x) 6K^{y)6xSJB{s) SA^{x) 

6K,iy)6xSJx{s) 5A^ix)6Axis) ^ ' ^ 



6K^{y)6xSJB{s) SA^{x)6B{s) 
When the relevant Fourier transform is taken and the result of (2.26b) is applied, the 
second term becomes 

^p^ f d^xd'^ye'P-''-'P-y{O\T{D,c{y)B{x)c{0)B{0))\O) 

= - / d^xe'P-''{O\T{B{x)B{0))\O) = 

where the ghost equation of motion, commutation relations, and (2.20a) have been applied. 
Thus after Fourier transforming (3.8a) we have 

F^Ap. -p. 0) = J'^xip. -p. 0) TxAp) (3.8c) 

where 

^.a(p,-P,0) = ^ J d^xd''ye'P-^-'P-y{0\TiDMy)Mx)mBm\0) (SM) 

is the full Green function. 

To one loop order, the contributions to this Green function are given by the diagrams of 
Fig. 2. The results of the calculation for the contribution to J-'^i^ from individual diagrams 
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of Fig. 2 are 



1T^ 



r,'p{ii 



9 'T'iiii 






2^i. 






The sum of these d 






g'N^ 



647r2p4 
' 647r2p4 

"647rV 

g'N, 



p^gt^v 



e 



-2e-|)+PMP.(-f + 6^+1 



P^fi-M^- ( ~ ~ 2 j + p^p^ ( - - - + 2 



647r2p4 
647r2p4 



PfiPi^ 
PiiPu 



^ 3 ; 

e e 



(3.9a) 



^ 3 ; 

T -T+4 
e e 



i 



agrams, including the tree level contribution is 



2J^M!^ = —. ^. o . (P £'m'- - PmP^) t + 2 + 2^ 



p4 647r2p4 
Note that we have introduced 



(3.96) 



1 1 



log 47r + 7 + log 



p. 



D = 4 + 2e . 



(3.9c) 



e e V '^ 

to be consistent withl^'^^, and included a factor of 2 to take care of the "crossed" term in 
(3.6d). 

It is important that this Green function is explicitly transverse beyond tree-level as 
required by conservation of the ghost current K^ = D^c, which follows from ordinary 
BRST symmetry. This is easily seen by considering the Green function J-'^,^. 



^q" I d^xd^ye'P-^+''i-y{O\T{A^{x)D,c{y)c{0)Bm\O) 



1 



d'^xe'P-''{0\T {A^{x)B{0))\O) ■ 



pp 



(3.9rf) 



where as usual, the ghost equations of motion and commutation relations have been ap- 
plied. Clearly the identity (3.9d) is satisfied by the tree level contribution p'^g'^/(2p2g2), so 
the higher-loop contributions to J^^u{Pi —p^ 0) must be transverse, serving as a consistency 
check on our calculation. 
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Returning to the identity (3.6d) and including terms up to one-loop we have the 
following result for Txf^{p) t^^l 



where 



U{p' 



9^ 
167r2 



13 A 1 97 1 ^2- 
Y~2J?~36~2^~T 



(3.10a) 



and some terms independent of the gauge parameter have been neglected. These neglected 
terms in (3.10a) come from quarks, and are irrelevant for the analysis of (3.6d). Combining 
(3.10a) with the results of (3.9) to one-loop order yields 



2F,Ap, -p. 0) = -{^^ if9,u - P.P.) [j + ^ + 2eJ 



(3.106) 



which confirms the QCD Nielsen identity (3.6d) for the gluon propagator. 

A formal solution for the gluon self-energy n(p2) can be obtained from (3.6) and from 



the transverse nature of F^^{p, —p, 0). Defining 

2F^^(p, -p,0) = {g^,.- 
then using (3.10) and (3.6) we have 



p2 



i.(l)(p2) 



(3.11a) 



^=[i + n(p^)]FW(p- 
^iog[i + n(p2)]=FW(p2) 

l + n(p2)=exp I F^^\p^)di 



(3.116) 



This formal solution relating the gluon self-energy to another Green function could be of 
interest in other contexts, but at present we do not see any immediate implications. How- 
ever, this result could provide a different approach for studying questions of confinement. 
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3.2 Gauge Dependence of the Quark Propagator 

The Nielsen identity for the quark propagator can be obtained directly from (3.5b) by 
functional differentiation with respect to i(^{x), i(^{y) and using quark number conservation. 

d d^r d^r s^T d^r sr 

+ 



5^T d^r d^r d^r " ^^'^^^ 

Si/j{y)dK^dx S'i/j{x)di/j 6-^{y)6-^ di/j{x)dK^dx 
As in Section 2, the first and second terms are zero because of Lorentz invariance with 
operator insertions at zero momentum. This leads to a result identical to (2.13) t^J. 

^S-\p) = S-\p) {F{p,-p,0) + F{p,-p,0)} . (3.13) 

Again, since S~^{p) is zero at the mass shell, and since F{p, —p, 0) has no single particle 
pole, we have the following Nielsen identity for the quark propagator. 



dC 



-ip) 



0. (3.14) 



Just as in Section 2, this result may be shown to imply that the quark mass shell is gauge 
independent in QCD for on-shell renormalisation schemes. We will return to a perturbative 
consideration of this idea in Section 4. 

3.3 Nielsen Identity for the Ghost Propagator 

To construct the Nielsen identity for the ghost propagator we need to differenti- 
ate (3.5a) with respect to c{x) and c(0). After applying ghost and fermion number con- 
servation the following identity is obtained: 

d d^r __ s'^T d^r _ s^T s^r 

d^Sc{x)6c{0) ~ ~ SxSKf,{w) Sc{x)6c{0)SA^{w) ~ SKf,{w)Sc{x) SxSc{0)6A^{w) 

6c{w)6c{Q) 5xSKc{w)5c{x) dxSc{w)dc{x)dc{0) 

The dependence on the variable w (which has an integration associated with it) has 
been made explicit to avoid possible confusion with c{x) . The first term on the right hand 
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side of (3.15) is zero since after Fourier transforming we have a zero momentum insertion 
with a vector operator. The third term is also zero since the one point function has the 
property B = dW/dJs = 0. 

To analyze the remaining terms we first recognize that 

[ d\x - y) e^P<^-y \ , ff^ , =^D-\p') (3.16) 

J 6c{x)6Kf,{y) p^ 

where D~^{p'^) is the inverse ghost propagator. After Fourier transforming (3.15) we have 
the following Nielsen identity for the ghost propagator. 

^,D-\p') = ^G^{p,-p,0)D-\p') + G^'\p,-p,0)D-\p') (3.17a) 
oE, p 

G'^'Hp,q,-p-q)^ I d'xd'ye^^-^+^^-y ^l^ (3.176) 

J 6Kc{x)6c{y)6x 

G^ip,q,-p-q)^ f d'xd'ye^^-^+^^-y /'j-^ (3.17c) 

J 5A^{x)5c{y)6x 

Now consider the Green function G^ to take into account the mixing between A^ and B. 



5xSc{y)5A^{x) J 6xSJc{w)J\{s) 6c{y) 6A^{x) 



d'^sd'^w 



6xSJc{w)SJb{s) Sc{y) SAf,{x) 



(3.18) 



SxSJc{w)Jx{s) Sc{y)6c{w) SA^{x)6Ax{s) 
+ I d'^sd'^w- 



,4 A 5'w^ <^'r 6^T 



dxSJc{w)dJB{s) Sc{y)Sc{w) SA^{x)SB{s) 

After Fourier transforming we find that the first term in (3.18) does not contribute to 
P'^G^{Pj —p, 0) since r^x{p) is transverse. Thus we have 

^G^ip,-p,0) = D-\p^)^ I d^xd'ye'P-^-'P-y{O\T{c{0)B{0)c{x)Biy))\O) (3.19) 

where (2.26b) has been used. 

Again, we can explicitly display the perturbative content of (3.17a) to one-loop by 
calculating the diagrams in Figures 3 & 4 and noting that the tree level contribution in 
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each case is zero. The results are 






P 



G(i)(p,-p,0) 
The result for the ghost propagator is ^^^' 



9^ 
647r2 

647r2 



1 

- -2 
e 



+ 2 



(3.20) 



^-i(p2)=p2[i + n(p2)] 

g^N, r 3 - ^ 



n(p^ 



167r2 



4e 



-1 



(3.21) 



Substituting (3.20) and (3.21) into (3.17a) verifies the Nielsen identity for the ghost prop- 
agator to one-loop. 

Similarly to the case for the gluon propagator, (3.17a) has a formal solution. Defining 



p^ 



G(2) (p, -p, 0) = G(i) (p, -p, 0) + ^G^ (p, -p, 0) 



p 



(3.22) 



and noting that there is no dependence on Dijp') in (3.22), we have the following formal 
solution to the Nielsen identity. 



D {p ) = p exp 



deG(2)(p,-p,0) 



(3.23) 



4. On-Shell Fermion Mass Renormalisation Constant 

In this section we wish to investigate the consequences of the Nielsen identities [see (2.14) 
and (3.14)] for the fermion propagator in a perturbative analysis of on-shell renormalisation 
constants. The key result from both QED and QCD is 






p2=M2 



0. 



(4.1) 



Consider the following definition of the fermion mass renormalisation constant (valid 
for both QED and QCD) in an on-shell scheme'^^'^^J 






(4.2) 
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where f = M is a zero of the inverse propagator. We will prove that Z^ (and hence M) 
must be gauge parameter independent to all orders in perturbation theory. 

To construct this direct link between the results of perturbation theory [^^~^^] and 
our formal proof of the gauge independence of the mass-shell we will show that the coeffi- 
cients Mi in the expansion of the mass renormalisation constant (we are using D — ^ — luj 
to be consistent with [9]) 



oo 



are gauge independent to all orders in perturbation theory. This has been explicitly ob- 
served to two loops in perturbation theory l^^"^^]. 

Now, the Feynman propagator is written, using the conventions of 1^^^, as 

Sr\v) = f -mo-E(p). (4.4) 

sw = g (^)" i^^' (f ) + (f - "'»>^- (f )) ■ (^-s) 

We must also have, in addition to the above Nielsen identity, that the inverse prop- 
agator itself must vanish on the mass shell (j? = M^). This is sufficient to define the 
coefficients Mj in the expansion of Zj^ in the following way. Working to the one-loop 
level, we will disregard anything of the order of a^ or higher. Thus we have the expression 
for Sp^{p) at one-loop to be 

S-.(p) = J, _ ™„ _ ^ (™„^, (^) + (^ _ „„)B, (^^)) . (4.6) 

We then substitute for mo from (4.3), obtaining, to order a^ 

s^\p) =f -m + m-m{i + ^Mi) 

This expression is then evaluated on the mass shell and set equal to zero: 

S-\p) ^_^MMi--^MAi(l)=0, (4.8) 
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thus we must have Mi = —Ai{l) in order for equation (4.8) to be satisfied. 

If one now considers the Nielsen identity (4.1), then one obtains to one loop the 
relation 

dS'^^ , . as f dAi f mf^\ , .dBi fmn\\ , ^ ^. 

-5f W = -^ ("o^ (^ j + (,. - ra,,)^ (^ j j . (4.9) 

Upon going on-shell and substituting for toq, we get 



dS^^ 



Thus we obtain the result that 



as ,.dAi 



F (p) ^__^M^(1)=0. (4.10) 



and thus that 



^(1)=0, (4.11) 



dMi 

^ 0, (4.12) 



dC 



ie., that Mi is gauge independent. 



It is this result that we wish to argue can be continued to all orders in perturbation 
theory. That this is possible can be seen by comparing the structure of the Nielsen identity 
above with that of the determining equation for Mi , before one goes on-shell. The structure 
is identical, with the exception of the leading term containing Mi. What remains when one 
goes on-shell, in (4.7) defines Mi, and what remains after going on-shell in (4.9) shows Mi 
to be gauge independent. 

Let us illustrate by continuing on to two loops. To the two-loop level, we have 



S± (™„^, (^) ^ ,y _ „„)B. [f 
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again we substitute for mo and keep all terms of order a^ or less. This gives 



S-\p) =p -M- f^MMi + -^MM2] ~^ImA 



P^ 



+ f # - M - ^^MMi) Si 



(4.14) 



^2^ V '\p^ 

Going on-shell, it is easy to see that the terms in (4.14) of order a^ are simply the one-loop 
result, and that the order a1 terms will define M2 to be: 

M2 = - (^2(1) + Ml (Ai(l) + 2A\{1) - Si(l))) . (4.15) 

On the other hand, we can see that application of the identity (4.1) to the two-loop 
expression (4.13) will result in 

dSp^ , , as f^^dAi {M^\ as _a /_ , (M 



2 



as 2M3 d / , ,, fM^W f , \ aSi fM^ 



a, _ d /_ ^ /M 
-M— MiBi 



(4.16) 



M-2^ d^ \ ' V p^ 

^BA2 fM^\ , ,,,dB2 fM^W 






and here we can see that when we go on shell the order as part reproduces equation (4.11) 
and the order a^ part will give 



F 



d^ 



(p) =~^{M^{A^{1) + 2A[{1)-B^{1)) + A2{1)) = 0, (4.17) 

which obviously guarantees that M2 is gauge independent. Note also that we made no use of 
one-loop results in obtaining equations (4.15) and (4.17), and thus the gauge independence 
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of M2 in no way depends upon the gauge independence of Mi, although the one- loop results 
will obviously make things simpler. 

One can see, then, that the relationship between the defining equation for M„ and 
the n-loop Nielsen identity (4.1), which shows the gauge independence of M^, will be 
maintained to all orders in perturbation theory. 

Having demonstrated that Z^ is gauge independent {dZm/d^ = 0) it then follows 
from (4.2) that dM/d^ = 0. 

dZm dM dZm _ Z^ dM 

^ ^ ^ (4.18) 



dM 



Z/iYi OZ/jji 



M dM 






This completes our perturbative analysis of the gauge independence of the mass shell and 
mass renormalisation constant in on-shell renormalisation schemes. 

5. Conclusions 

In this paper we have studied the Nielsen identities for the two-point functions of QED 
and QCD in the covariant formalism. It was demonstrated in the case of QED that the 
identities lead to results complementary to those of the usual Ward identities. As with 
the Ward identities, the Nielsen identities offer possibilities to check one's calculations, 
however, they also allow us to see where physical meaning may be found in apparently 
gauge dependent Green's functions. In particular it was proven that the electron pole 
mass, the photon polarization and the on-shell mass renormalisation constant Z^ are all 
independent of the gauge parameter. 

For QCD it was demonstrated that the quark pole mass and on-shell mass renormal- 
isation constant are gauge independent. This is a formal property and we are not sure 
as to its correct interpretation in relation to quark confinement. In particular the gener- 
ating function used above to describe QCD, does not take the Gribov ambiguityl^^^ into 
account. It has recently been shown that this ambiguity prevents the construction outside 
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of perturbation theory of a globally BRST invariant field with quark number one'^^l which 
could be identified with a physical quark. Perturbatively invariant quark fields do exist, 
however, and we stress the need to clarify the connection between the invariance of the 
pole mass, shown here using functional methods, and the perturbatively gauge invariant 
solutions with quark number one found in [15]. In this context we note that the gauge in- 
dependence of the pole mass has already been shownl^"^] to hold in perturbation theory up 
to two loops and appears to be true in the operator product expansion (OPE) of the quark 
propagatort^'^^], although care must be taken there since gauge-dependent condensates 
appear in the OPE of the gauge-dependent propagator [^^1. 

We have calculated the explicit one-loop content of the Nielsen identities obtained 
for the gluon and ghost propagators. Formal solutions to the ghost and gluon Nielsen 
identities were constructed for the propagators in terms independent Green functions. 
Since this gives a different view of the gluon propagator this solution may be of interest 
for studies of confinement. 
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Figure Captions 

Figure 1. Radiative corrections to tiie mixed propagator {0\T (B{x)Af^{y))\0). Ttie 
dotted line represents the field B. 

Figure 2. One-loop contributions to ttie Green function J-'f^^{p, —p, 0). Tfie representation 
of Afj_ and the ghosts is conventional, while the dotted line represents the field B. 

_ /J, 

Figure 3. One-loop contributions to the Green function ^G^{p, —p, 0). Diagram c leads 
to massless tadpoles and is thus zero in dimensional regularization. 

Figure 4. One-loop contributions to the Green function G^^^p^—p^O). Diagram b leads 
to massless tadpoles and is thus zero in dimensional regularization. 
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